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Abstract 

' As is well-known, if the scattering cross section of the Dark Matter 

(DM) particles off the nuclei inside the Sun is large enough, the Sun 
can trap DM particles. In principle, the annihilation of DM pair inside 
^ I the Sun can then lead to the detection of a relatively large flux of neu- 

trinos in the neutrino telescopes. If the annihilation directly produces 
a neutrino pair, the flux of neutrinos on Earth will be monochromatic. 
In this case, the oscillatory terms in the oscillation probability lead to 
a novel seasonal variation of detected events which is sensitive to the 
initial flavor composition. In this paper, we propose two models that 
predict such a detectable monochromatic neutrino flux from the DM 
annihilation. Model I, which is based on augmenting the type II see- 
saw mechanism, predicts a flavor composition for the monochromatic 
flux determined by (rn-jy)^^. In model II, the DM pair first annihi- 
lates to a pair of sterile neutrinos which oscillate into active neutrinos 
with a flavor composition determined by the flavor structure of the 
active-sterile neutrino mixing. 

1 Introduction 

Although reasonably large amount of evidence has been accumulating in fa- 
vor of Dark Matter (DM) as the explanation of the missing mass problem of 
the universe, the nature of DM particles is still unknown. Under the assump- 
tion that the DM particles are produced thermally in the early universe and 
there is only one kind of DM, the measured DM abundance in the universe 
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determines the self-annihilation cross section to be 1 ph. However, the DM 
mass or its annihilation products are still unknown. 

Arguably one of the most plausible classes of DM is Weekly Interacting 
Massive Particles (WIMP) with a mass, uidm, in the range 100 GeV — 
few TeV. This range of DM mass is interesting both for direct and indirect 
techniques of DM particle detection. Direct detection techniques are based 
on measuring the recoil energy in the scattering of DM particles off nuclei in 
a target. Indirect detection techniques are based on detecting the particles 
produced by DM annihilation in regions such as the Sun or the galaxy center 
where the DM concentration is relatively high. 

Within certain popular models such as MSSM, the DM pairs first anni- 
hilate to pairs such as bb, tt, ZZ and W^W~ . Neutrinos are then produced 
as secondaries through their decays. However, as discussed in [1], there are 
various ways to build a model in which DM particles dominantly annihilate 
to neutrinos. 

In this paper, we are interested in detecting the neutrino flux from the 
annihilation of DM particles that have been accumulated inside the Sun. 
These DM particles are non-relativistic. As a result, if neutrinos are directly 
produced by the DM pair annihilation, the spectrum will be monochromatic 
with El, — ruDM- However, if neutrinos are secondary products of the DM 
annihilation, their spectrum will be continuous with E^, < moM- In either 
case, if the DM mass, uidm, is larger than the detection threshold of the 
neutrino telescopes, a high energy neutrino flux pointing towards the Sun 
is expected at ICECUBE. Remember that the energy of ordinary neutrinos 
produced by proton fusion in the Sun center is too small to be detectable by 
neutrino telescopes with detection energy threshold of Eth few 10 MeV. 

The purpose of the present paper is to build models within which the 
method proposed in [2, 3] is effective. Such a model should have the fol- 
lowing properties: (1) The Dark matter mass is larger than the detection 
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energy threshold of the neutrino telescope: 0(50 GeV)-O(100 GeV). (2) The 
DM pair dominantly annihilates to a neutrino pair with a non-trivial flavor 
composition and a total annihilation cross section of 1 ph. (3) As it will 
be discussed in the next section, to obtain enough statistics, the scattering 
cross section of DM particles off nucleons should be greater than 10~^ ph. 
On the other hand, it should be smaller than 10~^ ph to evade the bounds 
from direct detection [4]. Other examples of such models can be found in 
Ref. [5]. 

In sect. 2, we briefly review the method proposed in [2, 3] and formulate 
the conditions a model has to satisfy to predict a sizable seasonal variation 
in the number of events from the DM annihilation in the Sun center. In sect. 
3, we propose Model I which embeds type II seesaw mechanism. Within this 
model, the flavor structure of a-(DM-|-DM — )■ I'o.i'/j) is determined by {m„)ai3- 
In sect. 4, we propose Model II within which the DM pair first annihilates to 
a pair of sterile neutrinos and then the sterile neutrinos oscillate into active 
neutrinos on the way to Earth. In sect. 5, we review our conclusions. 

2 A novel method to extract information about 
dark matter particles 

The DM particles propagating in the solar system have velocities about a 
few hundred km/sec. When these particles enter the Sun, they can lose their 
kinetic energy by scattering off the nuclei. They will then fall in the gravi- 
tational well of the Sun. As a result, during the Sun lifetime, the density of 
DM in the Sun has increased. The number density increases with the scat- 
tering cross section of these particles off nuclei. The trapped DM particles 
virilize and come to thermal equilibrium with the nuclei in the Sun center. 
Equating \Ekinetic\ = SkBTQ/2 with \Vgravity\ = 4:7TGNpmDMrl)M/^^ we find 
that the virilized DM particles are centered in the Sun within a volume of 
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radius tdm ~ {QkBTQ/SirGNPe'mDMY^'^- The DM particles inside the Sun 
are non-relativistic so the annihilation of a DM pair into a pair of on-shell 
particles will result in a monochromatic spectrum. In case that a neutrino 
pair is directly produced by the annihilation of the DM pair, the energy of 
each neutrino will be equal to the DM mass, moM- However, if the neutrinos 
are the decay products of unstable particles produced by the DM pair anni- 
hilation, their spectrum will be continuous. In the latter case because of the 
very large distance between the Sun and Earth, Lsun- Earth, the oscillatory 
terms in neutrino oscillation probability given by sm{Am'^Lsun-Earth/'^E^) 
will average to zero. However, as discussed in [2], in the former case where 
the spectrum is monochromatic, the oscillatory effects in the oscillation prob- 
ability are not averaged out and therefore lead to a seasonal variation in the 
number of events at ICECUBE as the Sun-Earth distance varies during a 
year because of the eccentricity of the Earth orbit. In [3], it was shown that 
studying this seasonal variation provides information on the flavor structure 
of DM -I- DM — >■ i'a+ p as well as on the value of moM (through the 
combination Am2i/mDM)- In particular, observing oscillatory behavior on 
top of the trivial variation of inverse of the square of the Earth Sun distance 
{i.e., L'^arth-Sun) ^leans the initial flux is monochromatic with a non-trivial 
non-democratic flavor composition {i.e., F^,^ : F^,^ : F^,^ 7^ 1 : 1 : 1). 

In [3], the details of this novel method to extract information on the 
properties of DM annihilation modes have been discussed. Of course, to study 
the time variation of the neutrino flux, sufficient statistics is also required. 
The statistics is determined by the rate at which the Sun traps DM particles 
and this rate in turn depends on the scattering cross section of the DM 
particles off the nuclei in the Sun. Scattering of DM particles off nucleons 
inside a nucleus is coherent so in the case that scattering is spin independent, 
the cross section of DM particles off a nucleus of Z protons and A—Z neutrons 
is proportional to \ZM.p + {A — Z)M.n\^ where M.p and M.n are respectively 
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the scattering amplitude of the DM particles off a single proton and neutron. 
(Notice that if the scattering was not coherent, the cross section would be 
given by Z|A^pp + {A — Z)|A^„p.) Because of the non-linear dependence 
on Z and A — Z, heavier nuclei can trap the DM more effectively than the 
same number of separate nucleons. Thus, although the majority of the mass 
of the Sun is composed of protons, in evaluating the capture rate via spin- 
independent scattering, the presence of the heavier nuclei inside the Sun has 
to be taken into account. It can be shown that for moM > 100 GeV, if 
the spin-independent cross section is larger than 10~^ pb and each DM pair 
produces a neutrino pair, a few hundred events can be registered each year 
in a detector such as ICECUBE [3, 6]. In this mass range, the direct bounds 
on the cross-section of the spin-independent scattering of DM particles off 
nucleons is ~ 10~^ pb. 

In [2], various possible sources of widening of the monochromatic spec- 
trum from DM-I-DM — > Ua + ^/s have been studied. The main source of 
widening turns out to be the thermal distribution of the velocities of the 
initial DM particles. This thermal distribution widens the monochromatic 
line to a narrow Gaussian with a width AE/E ~ 10~^(100 GeV/mDMY^'^ 
[2] . At these energies a fraction of active neutrinos undergo interaction with 
the nuclei inside the Sun while they cross the Sun. The charged current in- 
teractions of and produce charged leptons which are absorbed. The 
charged current interaction of produces a charged tau which in turn de- 
cays producing u ^ with a lower energy. The neutral current interaction of 
the neutrinos produces a neutrino of a lower energy. The outcome is that the 
scattering will reduce the height of the monochromatic line but will add a tail 
to the spectrum. Since the cross section of the neutral current interaction 
is finite in the forward scattering, the sharp line will remain sharp. At the 
sun surface the spectrum will be composed of a sharp line superimposed on 
the end of a continuous spectrum consisting of the scattered and regenerated 
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neutrinos. For m^M < 500 GeV, the energy and therefore cross section of the 
produced neutrinos are small which means a good fraction of the neutrino 
will remain unscattered resulting in a significant seasonal variation due to 
the oscillatory terms in the oscillation probability. However, for larger ttidm, 
the majority of the neutrinos are supposed to be scattered before leaving the 
Sun. This means in general for ttidm > 500 GeV, the method proposed in [3] 
is ineffective. In building Model I, we take this consideration into account. 
However, in the case of Model II, DM pair first annihilates into sterile neutri- 
nos that leave the Sun unscattered. These neutrinos then oscillate to active 
neutrinos. This means in the case of Model II, even if uidm > 500 GeV, the 
method might be effective. 

On the other hand, if itidm < a few x 10 GeV, the energy of neutrinos will 
be below the detection threshold of ICECUBE. Moreover, for lower values 
of the DM mass, the correlation between the direction of incoming neutrinos 
and the produced charged lepton will be lost. This in turn means using 
the directionality to reduce the background from the atmospheric neutrinos 
will become less efficient. Considering these facts, in model building we set 
ruDM > 100 GeV. 

3 MODEL I 

In this section, we build a model for neutrino mass and dark matter by slightly 
augmenting the type II seesaw mechanism. Through the standard type II 
seesaw mechanism, this model will lead to a Majorana mass for neutrinos. 
The DM candidate in this model can be either a complex scalar or a Dirac 
fermion that we add to the model. A Z2 symmetry guarantees its stability. 
As is well-known, within the type II seesaw mechanism, a SU{2) triplet 
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scalar exists with a nonzero hypercharge as follows 

(1) 



V2. 

The left-handed leptons couple to this triplet as follows 

(2) 

in which a and /3 arc flavor indices, i, j and A; are SU{2) indices and C is 
the 2x2 charge conjugation matrix: Cn = C22 = and Cu = —C21 — 1. By 
assigning lepton number equal to —2, the lepton number will be preserved 
by the coupling in Eq. (2). The pure scalar part of the potential involving 
only H and A is given by 

V = m^H^ -H+MlTrlA^A] + ^{H^- Hf + X^H^ A"^ AH +\;{rr[A^ A]H^ ■ H 



+^(Tr[AtA])V^|Tr[AA]r 



(3) 



+// [{Ay {H'^f + V2A- H+H"^ + A— . 

Notice that other forms of quartic couphng of H and A which preserves 
lepton number and electroweak symmetry can be rewritten as the combi- 
nations of the above terms. For example, Tr(AA''^[A''^, A]) = |Tr(AA)p — 
[Tr (At A)]2, Tr (At AAt A) = (Tr(AtA))2 - |Tr(AA)|V2 or H^AA^H = 
Tr{H^H)Tr{A'<A) - H^A^AH. If the quartic couplings are all positive the 
potential will be stable; i.e., as A and/or H ^ 00, V remains positive. 
In fact, a weaker condition guarantees stability: For example, as long as 
Aai > 0, the condition Aai + Aa2 > guarantees stability even if Aa2 is 
negative. 

The terms in the third line of Eq. (3) can be rewritten as tikHiA\-Hj. 
Since we have assigned lepton number of —2 to A, lepton number will be bro- 
ken only softly by /x. Along with H, A° also receives a tiny VEV proportional 
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to the lepton number violating /j, parameter 

Since fiv/m\ ^ 1, we expect (A°) ^ {H). The first term in Eq. (2) then 
gives a Majorana mass to neutrinos 



Although the coupling in Eq. (2) preserves the total lepton number, it 
violates lepton flavor and can therefore give a significant contribution to 
Lepton Flavor Violating (LEV) processes. Eor a comprehensive review see 
[7]. In particular, the bounds on — > eee sets the bound feefne/'iTT'A < 
1.2 X 10-5 TeV-^. Taking mA ~ 1 TeV and fap ~ 0.01 - 0.001, all the 
bounds will be satisfied. 

Eor mA ~ TeV and fap ~ few x 10"^, the values of // smaller than 10 keV 
result in small enough neutrino mass. The smallness of the /i parameter can 
be justified by 'tHooft criterion: In the limit that /i vanishes, the lepton 
number is preserved. 

Let us now discuss the DM sector. In the following, we discuss two kinds 
of candidates: (1) complex scalar; (2) Dirac Eermions. In both cases, we 
assign a lepton number equal to 1 to the DM candidate. We introduce a Z2 
symmetry under which only the DM candidate is odd. The Z2 symmetry 
stabihze the DM candidate. In order for the DM pair to annihilate to neutri- 
nos, we need to introduce another scalar, r], which mixes with A°. A similar 
model in the context of extra dimensions has been studied in [8]. Before 
discussing the DM particles, let us focus on 1]. We take to be a singlet 
complex scalar with lepton number opposite to that of A with the following 
interaction term 
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which after electroweak symmetry breaking leads to a mixing term: 



approximately equal to and m\ and mixing of X4V^/{2m\). 

A term such as rj'^H^ • if is a lepton number breaking term with a dimen- 
sionless coupling so we do not include it. We should protect rj from getting a 
large VEV; otherwise it will lead to a large (A°) and therefore large neutrino 
mass. A term of form 



leads to a mixing between rj and h given by ^'v/{m\ — mf^) and induces 
(77) — —^'v'^/m^. Notice that fi' breaks lepton number softly, so it should 
be also suppressed: /x' <C m^. The subsequent shift of A° will be given by 
A4(/x')^v^/m^. This small shift in (A°) does not change the situation. A 
small lepton number violating mass term of form rf can be also added to the 
Lagrangian but it has no serious impact on the discussion. 

As mentioned above, the DM can be either a complex scalar or a Dirac 
fermion. In both cases, the DM pair annihilates to a neutrino pair via a 
s-channel exchange of a scalar mass eigenstates that are mixtures of 77 and 
A°. Let us now discuss both possibihties one by one. 

• Complex scalar, $, as DM 

The general Z2 invariant and lepton number conserving Lagrangian 
involving $ can be written as 




Taking m\ ^ rri^ X^v^ /2, we shall have two mass eigenstates with masses 



(5) 





^($t . $)2 ^ ^^^^t . ^$t$ + X^^r]^r]^^ ■ $ + AA$Tr(AtA)*t • * . 
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The coupling in the last line mixes rj and A so leads to 

.(DM + DM = ^ _ (2„ _ . 

where = ""^1 + ^H>i>v^/'2. To obtain annihilation rate indicated by 

DM abundance in the standard thermal DM scenario, atot ~ 10~^^ cm^, 
the following relation should hold 

|™MW,)V 300 GeV'^i|(l^)\ (T, 

As mentioned earlier, a Z2 symmetry stabilizes $ against decay. If ($) 
is nonzero, the Z2 symmetry will be broken and various decay modes 
will become open for Vanishing ($) sets bounds on the parameters 
of the model. A necessary condition for vanishing ($) is 

ml^^ < 4(4A^$ + + X^){mjjj^ + m^) , (8) 

where is the quartic coupling of rj. If this bound is not satisfied, the 
minimum of the potential will lie at 

^ _ _ -3m^** ± [9m2^^ - 32(to| + m2)(A$ + A^ + 4A^$)]^/^ 

4(A$ + A^ + 4A^$) 

where + is for negative m^$$ and — is for positive mr^<^^. Eq. (7) 
combined with Eq. (8) imply 

— 2mDM\ < A4-\/A$ + A^ + 4A^$ GeV. 

This means a mild fine tuning between and 2mDM is required to 
obtain a tot ~ 10~^^ cm^. 

Lepton number violating terms such as 77$^$ and 77^$$ can be added 
to the Lagrangian but since their couplings should be much smaller 
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than m^$$, they cannot change the discussion. Moreover, once the 
lepton number is broken, a small mass term of form m|$$/2 can be 
also added. This means that there can be a splitting between the 
imaginary and real components of Notice that a real m^$$ coupling 
leads to the annihilation of a pair of real components together and a 
pair of imaginary components together. This is unlike the annihilation 
through a neutral gauge boson that takes place between the imaginary 
and real components of As a result, a small splitting between the 
real and imaginary components will not change the annihilation pro- 
cesses. However, the heavier component can decay into the lighter one 
and a pair of neutrinos via the mixing of A° and rj and via the fa/3 
couphng. If the decay takes place when the DM particles have become 
non-relativistic, the energy of neutrinos will be given by the mass dif- 
ference between lighter and heavier components of $ which is given by 
\rh%\/{2rnDM)- If this splitting is much smaller than 1 MeV (which is 
natural with taking m%/{2rnDM) ~ A* <^ 1 MeV), the energy of these 
neutrinos will be too small to destroy the products of nucleosynthesis 
even if the decay takes place at or after nucleosynthesis era. 

• Dirac fermion, ip, as DM 

The general Z2 invariant and lepton number conserving Lagrangian 
involving ip can be written as 

— 2 ^ '^s^ — 2 ^ 'moMtpRtpL + H.c. 

Again through a s-channel diagram, this Lagrangian leads to 

.(DM+DM = (,K-(2m,„)=lK-(2m„„)=lj 
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To obtain atot ~ 10 cm^, 



2 ^i^DM faf3 f 1 TcV 

500 GeV 0.01 V 



) 
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ml - {2mDMf\ 300 GeV' 
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Similarly to the case with complex scalar DM, to obtain the required 
value of (Ttot a fine tuning between and 2mDM is required. 

Lepton number violating mass terms of form mRip'^Cipii/2 and mLipJ^CipL 
can be added to the Lagrangian. These terms lead to a mass sphtting. 
Remember that we discussed the effects of decay and annihilation in the 
presence of mass sphtting for scalar DM. The same discussion apphes 
here, too. 

Notice that in this model the flavor structure of DM + DM — > i/qI/^ is 
given by Km,^)^^^. It is also possible to have a two component DM scenario 
with both ijj and A Z2 x Z2 symmetry can stabilize both of them. 

Through the coupling of the DM with 77 and the mixing of 77 with h, 
the DM can interact with nuclei but the scattering cross section will be 
suppressed by the lepton number violating parameter /i'; i.e., by the mixing 
of h and rj. To obtain significant scattering cross section off nuclei, we add a 
real scalar ^ as the portal to quark sector through mixing with the SM Higgs. 
The term that mixes H and ^ is 



The DM candidates also couple to ^. The couphng to the scalar DM, 



which leads to a mixing of 



tan2Q;;i^ = 
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leads to 

a(DM + N ^'DM + N) = a{^N ^N) 



fm^<s><^m^HH/ml'\^ f 2m GeVy [ GeVy [ fN 



""^ ^ ^ 0.1 J 1^ ruDM ) y ruh ) 
where N collectively denotes nucleons {i.e., n and p) and /ibmn — 1 GeV is 
the reduced mass of the ^-N system [9] . 

The same coupling and mixing lead to $ + ^ /i* ^ / + /, + 
W~, Z + Z with (av) equal to 

2 2 2 

(2mDMr(/i ^ final states)) |„i^-H.2moM 



4mV(4ml)M - "^|)^(4m|,^ - m2)2- 

(10) 

r(/i — >■ final states)) versus can be found in [10]. Taking m^$$m^ijij/|4m|,- 
m|| ~ 0.1 and uidm < 300 GeV, (atotv) will be smaller than 10"^^ cm^. How- 
ever, for niDM ~ 200 — 400 GeV, the annihilation mode via s-channel h and 
^ exchange can be significant along with $ + $ — > i/ + i/. As discussed in 
[3] , in this case still the seasonal variation of muon track events at neutrino 
telescopes can be significant. 

Similarly, the coupling to the fermionic DM, ijj, 

y^^^i^^ + H.c. 

results in 

a(DM + AT ^ DM + TV) = a{^N ^ V^TV) = 2^ (^ ^^^^f j /fDMiv< 



4 ~ 

mt 



10"^ ph 



Y^^iu.nii x3TeV \ / 120 GeV V //jy 
m| j \ ruh ) VO-3 



(11) 

2 
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where fiBUN — 1 GeV is the reduced mass of the ip-N system [9]. In both 
cases, 0.14 < /„, /p < 0.66. Thus, within the favored range of parameters 
{iTLoM ~ 200 GeV), the scattering cross section will be of order of 10~^ pb 
which is high enough to lead to a few hundred neutrino events (or even more 
for larger uidm) at ICECUBE each year. We can have a detectably large 
monochromatic neutrino flux from DM annihilation inside the Sun with a 
non-democratic flavor composition determined by {'m^)ai3- The Y^^ coupling 
also leads to + t/j ^ f + f, W+ + W-,Z + Z with (av) equal to 

{2mDMT{h ^ final states)) U,^2moM 77-2 m'^\2(Arr,2 :;W^''rel 

(12) 

where Vrei is the relative velocity of the DM pair. Notice that in Eq. (10) 
for the scalar DM case, such a factor of v^^i does not appear. At the decou- 
pling era, Vrei ~ l/-\/20 so for V^mDMiTi^HHY^^/\4:ni'j}j^ — m|| ~ 0.1, the 
annihilation to the Higgs decay products in the early universe can be only 
subdominant: {a{'ip + 'ip ^*,h* a,nyth.mg)v) / (atotv) < 0.01. For the DM 
particles trapped inside the sun, 

2 3A;bT -^Q-gSOO GeV keT 



u 1 J- vj 

TTlDM TTlDM 100 eV 

so the annihilation to the Higgs decay products can be safely neglected. Thus, 
for the purpose of this paper, the vf^i dependence is favored. Had we taken 
the interaction of ^ with ^p to be of form i^ip^y^i/j instead of ^i/jiIj, such factor 
of v^^i would not have appeared in Eq. (12). 

Let us now discuss the flavor structure of the neutrinos produced by the 
DM annihilation. The amplitude of DM + DM — )■ z/q, + z/^, Aia/B, is propor- 
tional to fai3 which is in turn proportional to {171,^)0(5 ■ This means in the mass 
basis where (mj,) is diagonal, the amplitude Ai is also diagonal. Let us denote 
the neutrino mass eigenstates in vacuum by \i) = |1), |2) and |3). The neu- 
trino production is in the form |li), |22) and |33) (but not for example of form 



14 



|12)). The production rate of is given by \fii\^ oc |mjp. That is the den- 
sity matrix in the mass basis is proportional to Diag(m^, ""^D- As a result, 
for the quasi-degenerate neutrino mass scheme with |mip ~ |wi2p — 1^3 
the neutrino production will be democratic and as discussed in detail in Ref. 
[3], the seasonal variation will vanish. Let us now evaluate the seasonal 
variation for the general neutrino mass scheme. Due to the matter effects 
inside the Sun, a pure mass eigenstate while crossing the Sun converts to 
a combination of mass eigenstates as follows: |?'; surface) — ^jCiij\j) and 
|i; surface) = Sj%|j) where due to the matter effects a^j may differ from 
ciij. ttij and Qij are unitary matrices: J^j'^ij'^lj — ^'^'^ '^j'^ij'^lj — ^ik- 
After traversing the distance between the Sun and Earth, L, \j) and \j) will 
pick up a phase of g"*'"?^/^^. Neglecting the average of the phase e*^'"3i^/2^^ 

we can write (P(i/, ^ v^)) = E,- l^^wf l%f + SafJfea.sC/^it/^se'^™^^''/'''] 
(see Ref. [3]). Considering the fact that the production of \i) is given by m?, 
the seasonal variation due to Aml2L/2E ~ 1 can be evaluated as 

variation = — 5- — i^.^, .„ . 

Since Am|i <c; Am|i, we can take ml ~ and can therefore write 

variation = 2 1 \2\tt \2 (13) 

where we have used the unitarity of a^-. If ^13 is exactly zero, despite the 
matter effects, |3) will not change which means a^i = 032 = 0, so the vari- 
ation will vanish. Remembering that the matter density of the Sun falls as 

g-r/(o.iii.un) find that for sin ^13 ~ [(0.1i?sun)V;|cenier]"^ ~ 0.001, thc val- 
ues of 031 and 032 can in general be of order of 1 which means the variation 
can be sizeable (larger than 10 %) and can be measured by a few hundred 
events. Eq. (13) also confirms that for quasi-degenerate mass scheme with 
ml ^ |Am|i|, the variation is suppressed even for 03^032 ~ 1/2. 
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Let us now briefly discuss the impact of tliis model for the collider searches. 
Like the case of type 11 seesaw mechanism, we expect a triplet scalar which 
is produced by electroweak interactions. Like the type 11 seesaw mechanism, 
the A"*" and A"*""*" components decay into lepton pairs via the fap coupling. 
For a review of the possible decay modes of A, see [11]. However, there will 
be a new decay mode for A°, which unlike A'' — )■ UaJ^/s, is visible. The new 
decay mode is A° — >■ rjh. If rj is heavier than 2m dm, it can decay into a 
DM pair which appears as missing energy. Since in this model H mixes with 
^, the signature of Higgs can be in principle different. In fact, any process 
studied for the SM {e.g., gg h and h — )■ rf, 77) should be reconsidered 
for two SM-like Higgs fields {i.e., two mass eigenstates composed of h and ^) 
with rates suppressed by cos^ and sin^ a^^. The production of the second 
Higgs-like scalar will be suppressed by a factor of sin^ a^^. For sin ~ 0.1, 
which is the favored range by bounds on a(DM + N ^ DM + N), the pro- 
duction of the second Higgs can be neglected. In this case, the Higgs sector 
will be similar to what we had within the Standard Model. 

4 MODEL II 

In this section, we introduce another model within which DM is composed of 
Dirac fields, ■0, which are protected against decay again with a Z2 symmetry. 
We introduce a new C/(l)' symmetry with gauge boson Z' under which all 
the SM particles are invariant. In this model, us with masses close to those 
of active neutrinos also exist such that oscillation can take place between the 
sterile and active neutrinos. 

The fields that are added to this model are the following: 

• A Dirac field, ip = {ipL V'r), which plays the role of the DM candidate; 

• One (or more) left-handed sterile neutrino, 1/3] 
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• Right-handed neutrinos ur/s which are singlets under both electroweak 
symmetry and U{iy. These neutrinos are included to give Dirac masses 
to the rest of neutrinos. 

• U{iy gauge boson, Z'; 

• A complex scalar field H' which is an electroweak singlet but under 
U (1)' has a charge equal to that of ipL- H' receives a Vacuum Expecta- 
tion Value (VEV) of v' / \/2 which breaks the U{1)' symmetry and gives 
a mass of m,z' = g'v' to Z' . After electroweak and symmetry 
breaking, the H' also mixes with H and acts as a messenger to interact 
with nuclei. 

Via a s-channel Z' exchange, the DM particles annihilate into sterile 
neutrinos which in turn oscillate into active neutrinos with a non-trivial flavor 
composition given by the flavor structure of the mixing of the sterile neutrinos 
with the active neutrinos. 

We assume that the right-handed is neutral under U{X)' and only left- 
handed couples to Z': 

With this charge assignment there is no U{1)' anomaly. If we assigned the 
same C/(l)' charge to ■0^, we must add additional chiral fermions to cancel 
the U{1)' anomaly. 

Through this coupling, the DM pair annihilates into a sterile neutrino 
pair with the annihilation cross section 

where the mass of us which is taken to be of order of the active neutrino 
masses is neglected {i.e., rriyg <^ moM)- To have {atotv) = 10~^^ cm^ {i.e., 
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the value suggested by the DM abundance in thermal DM production sce- 
nario), we find 

m|, = 4m^Af - (^) "^-DM X (100 GeV) . 

For example, for g' — 0.16 and vriDM — 200 GeV, the desired abundance of 
DM can be achieved with mz' = 375 GeV. The Z' boson in this model does 

not couple to quarks or Icptons and there is no mixing between Z' and the 
SM gauge bosons. As a result, none of the bounds from the new gauge boson 
searches at colliders applies in this case. 

After breaking of U{1)' and the electroweak symmetry, the neutrinos and 
ijj receive Dirac mass which comes from the following Yukawa couplings: 

>C = fa0^R^LiHhij + fR^PR^H'us + Y^i^R'^LH''^ + H.c. . 

We assumed that only is odd under the Z2 symmetry so the Z2 symmetry 
prevents a coupling of form pRpipLiH'^ or a mass term of form u'^ijCipR. The 
masses of ■^l and ipR are the same and equal to 

v' 

After H and H' obtain vacuum expectation values, vg receives a Dirac mass 
and mixes with active neutrinos. For illustrative purposes, let us consider 
one sterile neutrino, one active flavor and two right-handed neutrinos. The 
masses will be given by 



flaV /rW 






hocV fR2v\ 







(14) 



The mixing between sterile and active neutrinos is given by the following 
formula 

tan 29 — 



ilL + fL)v' - (/ll + fR2> 



.12 
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The Lagrangian of the model includes a term as follows 



In the unitary gauge H ^ {0 {v + h)/V2) and H' = {v' + h')/V2: 

1 [h h'] L< , ^^'^f] . (15) 

2 ^ ^ \Yh'hVv' ml, J \h'\ ^ ' 

The mixing between h and h' is given by 

2Yh:hVv' 



tan 2ahh' 



ml - ml 



This mixing leads to DM scattering off the nuclei in the target of the direct 
DM search experiments with a cross section 

^(DM + AT ^ DM + TV) = f ^^$^1 ' /^dmat^ = 

Y^Yh>hv'lmn'V /200 GeVV /120 GeVy //jv 



\ 0.1 J \ mh' J V mh ) 
As expected, this cross section is the same as the one in Eq. (11) provided 
that we replace m^ — >■ m^/, m^HH Yh'ht^' ^-nd Y^^ — >■ Y^/\/2. As men- 
tioned before, 0.14 < /„, fp < 0.66. The cross section is large enough 
to lead to a significant capture rate and subsequently to a large neutrino 
flux detectable at ICECUBE. Again replacing m^ — )■ m/i', m^uH Yh'hv' 
and Y^^ Y^/y/2 in Eq. (12), the formula for a{DM + DM -> C,h* 
W~W^, ZZ,tt) can be obtained. Similarly to the case of Eq. (12), we find 
that for Yhh'Y^niDMv' /\4:m%j^ — m\i\ < 0.1, this annihilation mode can be 
safely neglected. 

Among the new particles that arc added to this model, only h' has a 
significant coupling to the SM. As a result, only h' might appear in the 
collider searches through mixing with the Higgs field. The same discussion 
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in the case of Model I described in the previous section apphes here, too: If 
sin a^hi < 0.1, the production of the new Higgs at coUider will be suppressed 
by sin^ a hh' < 0.01 relative to the ordinary SM Higgs. 

Let us now address the oscillation of sterile neutrinos to active neutrinos 
and the flavor composition of the neutrino flux from DM annihilation when 
they reach the detector. This problem has recently been addressed in [12]. 
In Eq. (14), for simplicity we have assumed that us mixes only with one 
active flavor. In general, z/5 can mix with more than one flavor. In fact, 
explaining the LSND and MiniBooNE results in the context of the present 
scenario suggests that the sterile neutrino simultaneously mixes with z/g and 
Vfj,. Recent fits to the short baseline neutrino data can be found in [13]. 
Some hints for the oscillation of atmospheric to us has also been found in 
the ICECUBE data [14]. The scenario proposed in this section can include 
two sterile neutrinos mixing with different flavors as required by the 3 + 2 
scenario suitable for explaining the LSND and MiniBooNE results. Some 
hints for sterile neutrinos have also been found by studying the solar neutrino 
oscillation data [15]. Studying the general case is beyond the scope of the 
present paper. We focus on the particular case that the mixing matrix in the 
(fe ly^ fr ^s) basis is of the following form 

=034-(t/pMiV5®llxl) 

where 

'1 ■ 

_ 1 
cosOu sine34 
—sin ^34 cos ^34 

{UpMNS © lixi) is the external sum of the standard three by three mixing 
matrix of neutrinos Upmns with lixi which means its flrst 3x3 block is 
UpMNS-i its 44 element is equal to one and the rest of its elements vanish. 



(16) 
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Notice that, without matter effects, P(i/e— ^-i^e), P{i'ix^i'fj,) and P(i/^^i/e) 
as well as P(^Ue^i^,i) are the same as what we expected without the ster- 
ile neutrinos. As a result, the bounds from reactor searches as well as the 
measurement of P(i'^^— ^i/j) and P(i'^^— ^i/^^) by atmospheric and long baseline 
experiment and various other experiments do not apply here. This mixing 
affects P{uh^Ut) which is poorly constrained by observation. The strongest 
constraint comes from the measurement of the total neutral current inter- 
action of the beam in the MINOS experiment [16]. If this constraint is 
saturated, ^34 can be still relatively large leading to a significant probability 
of sterile to active oscillation. 

A part of the us flux will oscillate into active neutrinos on the way to the 
Earth which can be detected at ICECUBE or its deepcore (provided that 
the detection threshold is lowered below uidm)- The number of muon track 
events in a time interval is given by {P{j^s ~^ ^fi)) which is the oscillation 
probability averaged over the time interval. As discussed in [3], due to the 
averaging effects the oscillatory terms given by Am|]^ are subdominant how- 
ever we should keep the oscillatory terms given by Am|]^. A simple numerical 
calculation shows that the matter effects on the oscillation probabihty cannot 
be neglected. Adopting the formalism in [3] and generalizing it to a four- 
neutrino scheme, we can expand the evolved 1/5 at the Sun surface in terms 
of the mass eigenstates as follows 

li/5; surface) = as\\l) + a52|2) a53|3) a54|4) . 

By a numerical calculation taking into account matter effects, asi can be 
found in terms of the U^^^ elements and Am34. Averaging out the oscillatory 
terms given by Amg^, we can then write 

{P{vs ^ i^^^)) = Ki\^\asi\^ + |f/^2|>52r + \U^3\^\as3\^+ 

ocnrT-T- 7-7-* * / . A77l2i-£'Sun-Eaj:th \i 

2U[U^iU^^asias2eM^ t^b )] ' 
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where E„ = itldm- As discussed in [2], by studying the periodicity of muon 
track events, the value of Amf2/m£,M can be derived. The resuhs of a simple 
numerical calculation shows that if the bound from [16] on 634 is saturated, 
P{i^s ~^ can reach as large as 0.2. This means for o"(DM + iV — )■ DM + 
N) ~ 10~^ cm^, up to a few hundred muon-track events can be registered 
each year at ICECUBE so the statistics will be enough to study a seasonal 
variation [3]. In ref. [3], we have found that for three- neutrino scheme in 
the case that DM + DM — t- z/^ + z/^ or Ur + z/^, the variation will be small. 
Obviously this consideration does not apply to the present model with four 
neutrinos and DM + DM — >■ 1/5 + 1/5. 

Notice that {P{vs — >■ ^n)) is not sensitive to Am|4 because we have taken 
U,j.A — 0. The total oscillation probability to active neutrinos, (1 — P{vs — > 
^s)), is given by 

Wsi^lasi? + \Us2?\as2\'' + Wsz^las?,^' + \UsA?\asA?+ 
21k[UsiUs2asias2 exp(z — )] + 

o^rrr TT* * / .A^43 -^Sun-Earth 

2SR[?753f/54%3a54exp(z — )] . 

The number of the cascade- like events is given by (JArc(l ~ P{^s ^s)) + 
ccc(l ~ P{^s ^s) ~ P{^s ^m))- Thus, by studying the variation of 
cascade-hke events, it will be in principle possible to derive /^rn^^. The 
possibility can be realized only if the following three conditions are ful- 
filled: (i) The energy threshold for cascade detection is below ttidm', (ii) The 
statistics is high enough; (iii) The oscillation length Lose — 47rmDM/Am|4 
is of the order of the seasonal variation of the Sun-earth distance: i.e., 
Lose ~ ■A-^Sun-Earth ~ 5 million km, which means 
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It is noteworthy that for itidm ~ 200 GeV, this means Amf2 ~ Am|4. UnUke 
the case of /i-track events, the directionaUty cannot be used to reduce the 
background of shower-hke events from atmospheric neutrinos. At ICECUBE, 
we expect a maximum of a few hundred cascade-like neutrino events from 
the DM annihilation inside the Sun per year which is within the statistical 
fluctuation of the background from the atmospheric neutrinos. This means 
that after 10 years of data taking by ICECUBE, the confidence level of the 
discovery of such cascade-like signal will be about 3-4 a. To make a serious 
measurement, a detector ten times larger than ICECUBE with threshold 
below m£)M is required. After about 2-3 years of data taking by such a 
detector, the statistics can be enough to claim discovery and after ten years, 
enough data can be collected to study variation and to extract information 
from the variation. 

Let us now discuss the possibility of constraining this scenario by vari- 
ous other oscillation observations. As discussed above, P(Pe — > Pe) relevant 
for the reactor neutrino data is not affected by the presence of the sterile 
neutrino. However, because of the matter effects inside the Sun, the sur- 
vival probability of the solar Ug will be affected by the presence of the sterile 
neutrino. The effect in general is rather complicated. For our case with 
|Am|i| ~ |Am|i| » Am2i, the simplified formalism in [17] can be appfied. 
Using the formalism in [17], we find that the deviation from SM neutrino 
oscillation without sterile neutrino is suppressed by (A^„/2A^e) cos2^i2c|3S34 
where Nn and N^. are respectively the neutron and electron number densities 
in the Sun. Inserting the numerical values, we find that the deviation is of 
order of 0.02 sin ^34 and is therefore very small even for sin ^34 1. 

Measuring P(i'J— ^i/t-^) at long basefine experiment OPERA or by study- 
ing the atmospheric neutrino data at Super-Kamiokande constrains ^34 as 
in this scenario F^V^^i^r) is suppressed by a factor of cos^ ^34 relative to 
the standard three-neutrino scheme. This scenario can be also probed by 
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improving the measurement of the rate of neutral current interaction events 
at long baseline experiments such as MINOS. 

5 Summary 

We have introduced two models within which the DM pair dominantly an- 
nihilates into a neutrino pair with a non-trivial flavor composition. As a 
result, the non-relativistic DM pair annihilation will lead to a sharp line in 
the neutrino energy spectrum with : : ^1:1:1. As shown in 

Ue Vt 

[2, 3], despite the very large distance between the Sun and Earth, the oscil- 
lation probability does not average to zero, leading to a seasonal variation as 
the distance between the Sun and Earth varies during a year. The variation 
contains information on vh^im and the initial flavor composition. 

In these models, the DM particles are stabilized by a Z2 symmetry under 
which only the DM candidates are odd. In both models, the DM particles 
interact with nuclei via a new scalar which mixes with the SM Higgs so the 
interaction is spin-independent. The parameters are chosen to yield a scat- 
tering cross section in the range 10~^ — 10~^ pb. Thus, while the present 
bounds from direct searches are satisfied, still significant statistics (a few 
hundred events per year) are expected at ICECUBE for indirect DM detec- 
tion. This means by a slight improvement in both direct and indirect DM 
searches, these models can be tested. The same couphng and mixing can also 
lead to DM DM j ,Z ^ Z, W . However, these annihilation 

modes will be subdominant. When DM is composed of Dirac fermions, these 
modes, being P-wave effects, are further suppressed by v^^i and for the case 
of the DM pair trapped inside the Sun can be safely neglected. 

Model I embeds type II seesaw mechanism so neutrinos are Majorana 
particles. Within this model, the DM candidates can be either complex 
scalars or Dirac fermions with lepton number equal to 1. It is also possible 



24 



that both the complex scalars and the Dirac fermions contribute to the DM 
in the universe. The DM pair annihilates into z/qZ^/j with a flavor composition 
determined by {171^)^13. The prediction of the model for LFV rare decays as 
well as the accelerator searches is similar to the predictions of type II seesaw 
mechanism, except that here the neutral component of the triplet, A°, can 
have a new decay mode, A° h + missing energy. 

In Model II, the DM is composed of Dirac fermions which via the exchange 
of a new Z' gauge boson annihilate into a pair of sterile neutrinos. Since Z' 
does not couple to quarks or ordinary leptons, it cannot be produced in the 
lepton or hadron colliders. In this model, neutrinos are Dirac particles so we 
expect a null result in searches for neutrinoless double beta decay. 

Since the annihilation products in the Sun center are sterile, they do not 
scatter off the nuclei present inside the Sun so, to the first approximation, the 
height of the sharp line in the spectrum is not reduced by scattering. On their 
way to Earth, the sterile neutrinos oscillate into active neutrinos. The flavor 
composition of the flux on Earth is given by the mixing parameters of the 
sterile neutrinos with the active neutrinos. The number of the muon-track 
events is given by P(i^5 — >■ i^^) which in turn is given by the active-sterile 
mixing. As discussed in the text, the main constraint on the mixing comes 
from the MINOS measurement of the total neutral current interaction of the 
beam at the detector. Taking into account this bound, the statistics can be 
still high enough to employ the method introduced in [3] . By improving this 
bound, the model will be further constrained. While the variation of the 
muon-track events at neutrino telescope will be sensitive to /Sm^^l'^HM-, the 
variation of the cascade events at the neutrino telescopes will be sensitive to 
both Am34/mDM and iS.m^^l'^DM- 
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